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*“■  A new,  potentially  disruptive,  phenomenon  in  target  pellets  irradiated  by  laser  or  particle  beams 
is  described.  Initially  weak  standing  sound  waves  (vibrations)  in  an  imploding  spherical  shell  can 
become  comparable  in  amplitude  with  the  shell  thickness  if  the  effective  adiabatic  index  satisfies 
" y < 5/3.  Detailed  calculations  with  a soluble  self-consistent  model  show  the  relative  amplification 
is  stongest  for  long-wavelength  modes  in  thin  shells  with  small  y,  and  is  insensitive  to  the  density 
profile.  I 
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Linear  Stability  of  Self-Similar  Flow: 

3.  Compressional  Waves  in  Imploding  Spherical  Shells 


When  a pellet  or  shell  of  a medium  containing  standing  sound  waves  undergoes  compres- 
sion, the  external  force  does  work  in  compressing  the  waves  as  well  as  the  medium  itself.  The 
wave  amplitude  tends  to  decrease,  but  the  characteristic  thickness  d of  the  compressed  medium 
may  decrease  faster.  When  this  happens  a/d,  the  amplitude  of  the  oscillations  relative  to  the 
size  of  the  system,  increases.  If  aid  becomes  of  order  unity,  the  perturbations  disrupt  the  basic 


We  can  determine  the  condition  for  such  "relative  instability"  using  the  following  simple 
argument,  applicable  to  systems  having  spherical  symmetry.  For  motions  in  which  the  external 
pressure  pit)  does  not  increase  too  abruptly,  the  action  aikc  is  an  approximate  adiabatic  invari- 
ant, where  k is  the  wave  number  and  c is  the  speed  of  sound.  For  adiabatic  compression, 

y~* 

c = (yp/p)'  1 — p ly  . The  characteristic  scale  size  d is  related  to  p by  d — p 1 3 — p 1 and 
to  k by  kd  — const.  Combining  these  expressions,  we  find 


a /d  — (dc)  1/2  — p 4y  . (1) 

Thus  if  y is  greater  than  5/3,  the  "geometric"  ratio  of  specific  heats,  the  relative  amplitude  of 

the  perturbations  decreases.  If  y < 5/3,  the  relative  amplitude  grows.  Similar  reasoning,  eg. 

in  slab  geometry,  leads  to  the  relative  growth  condition  y < 3.  It  should  be  emphasized  that 

the  mechanism  described  here  does  not  depend  on  the  magnitude  of  the  effective  acceleration  — — 

Section  1 

and  is  unrelated  to  the  Rayleigh-Taylor  instability.  0°^  Section 
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If  y — 1 (the  worst  ease),  ahl  — p'  b — pl/6.  Though  weak,  this  dependence  is  enough 
to  bring  about  a state  of  catastrophic  vibration  in  a target  pellet  imploding  under  the  action  of  a 
laser  or  particle  beam,  where  compression  by  as  much  as  106  can  take  place.  The  crucial 
question  is  the  appropriate  form  for  the  equation  of  state.  Thermal  conduction,  incomplete 
molecular  dissociation,  and  the  energy  required  for  ionization  tend  to  reduce  the  effective  adia- 
batic index  y below  its  geometric  value;  condensation  at  high  densities  tends  to  raise  it. 

The  aim  of  this  letter  is  to  illustrate  the  above  ideas  by  some  analytic  calculations  of  the 
amplification  of  acoustic  modes  in  imploding  spherical  shells.  For  this  purpose  a self-consistent 
soluble  ideal  fiuid  model  previously  used  by  Kidder  [1]  is  employed.  The  equation  of  state  is 
taken  to  be  a polytropic  adiabatic  law  with  arbitrary  y.  The  results  indicate  that  amplification  is 
determined  by  the  value  of  y,  not  the  form  of  the  density  profile. 

Consider  an  initially  stationary  spherical  shell  with  density  p(r)  in  the  region  defined  by 
#■_</■<  r+.  A pressure  loading  p+(t)  applied  unformly  at  r = r+  causes  it  to  begin  implod- 
ing radially  with  velocity  u(r).  Assuming  the  shell  material  to  be  an  ideal  fluid,  the  motion  is 
described  in  Lagrangian  variables  by  the  equations 


i>  + f>R~2  (R2u)  =0; 

(2) 

(3) 

(pp~>)'  *•  0, 

(4) 

where  R(r,  i)  is  the  position  at  time  t of  a fluid  element  whose  initial  position  was  r,  and  time 
derivatives  are  denoted  by  dots. 

The  class  of  solutions  of  Eqs.  (2-4)  which  describe  what  is  known  as  homogeneous  self- 
similar motion  results  from  assuming  the  existence  of  a function  f(t)  such  that  for  any  fiuid 
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with  ./'(0)  = 1,  /( 0)  = 0.  From  Eq.  (5)  it  follows  that  u(r,  I ) — r/,  while  the  solution  of  Eq. 
(2)  is 


p(r,  f)  — p0 (')/  3-  (6) 

Hence  from  Eq.  (4),  the  pressure  must  be  of  the  form 

p(r,  t)  =-  s(r)f>y  = s(r)ptff~3y  (7) 

For  self-consistency,  the  driving  pressure  must  balance  the  pressure  at  r = r+  found  from  Eq. 

(7),  i.e.,  />,.(/)  = p(r+,  i ) — s(r  +)pg(r  +)f  ~iy.  Substituting  Eqs.  (5-7)  in  the  force  law  (3),  we 

obtain 


- ff" 


+i 


.-2 


1 dp0 
p0r  dr 

where  a = 3(y  — 1)  and  r is  a separation  constant. 


1 ji 

t>or  dr 


(spg) 


(8) 


The  time-dependent  part  of  Eq.  (8)  can  be  integrated  once  to  give 


for  y > 1,  and 


r / = - 


— (/“'*  - 1) 
a 


1/2 


(9) 


rf  = - (2  In/'1) 1/2  (9) 

for  y - 1.  Equations  (9)  and  (9  ) distinguish  the  motions  called  Type  II  by  Sedov  (2),  who  first 

investigated  them.  Both  equations  evidently  give  rise  to  functions  / which  vanish  at  some  finite 

time  t0,  when  the  system  becomes  singular  (implodes  completely). 

The  spatial  part  of  Eq.  (8)  yields  a solution  for  each  non-negative  choice  of  s or  of  p0.  We 
select  two  cases  for  investigation:  (i)  Isentropic  motion,  s =*  s =»  const;  (ii)  Uniform  density, 
p0  = p = const.  Isothermal  (y  = 1)  motion  is  distinguished  as  subcases  (i)'  and  (ii)',  respec- 
tively. In  case  (i)  we  find 


Po  = P 


p(y  - l)r; 
2pyr2 


l/(y-D 


(10) 
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Here  p and  p = spy  are  constants.  In  case  (i)' 

P0  = p exp[prJ/(2pr2)].  (10) 

For  case  (ii)  we  find,  independently  of  y, 

s = s[pr2/(2pr2)  - lj.  (11) 

In  order  to  restrict  the  problem  somewhat,  we  choose  ri  = 2ypr2/(y  - l)p  in  case  (i), 

making  p<>(r_)  = 0,  and  rl  = 2 pr2/p  in  case  (ii),  making  s(r_)  = 0,  so  that  in  both  cases 
P-  * p(r-)  — 0.  Then  we  express  r in  units  of  r_.  For  case  (i)\  p(r)  nowhere  vanishes,  so  we 
take  r_  = 0 and  normalize  r with  respect  to  (pr2/p)1/2.  Time  is  conveniently  expressed  in  units 
of  t and  density  in  units  of  p. 

We  begin  by  investigating  linear  stability  with  respect  to  compressible  perturbations  for 
case  (i).  It  is  convenient  to  employ  the  formulation  developed  by  Book  and  Bernstein  [3,4]. 
The  fluid  element  whose  unperturbed  motion  is  described  by  R(r,  i ) is  displaced  to 
R(r.  r)  +Hr,  i ).  Expressing  the  perturbed  density  and  pressure  in  terms  of  the  Lagrangian 
displacement  ^and  the  unperturbed  quantities,  we  obtain 

r*i*vfe(v  Vln',»l|  (l2) 

-i  v v ln  Po  ' 

The  free  surface  boundary  condition  associated  with  Eq.  (12)  is  just  V • £ =»  0 at  r = r±. 

In  Eq.  (12)  the  space  and  time  dependence  are  completely  separated.  An  equation  for 
a = ^ X is  derived  by  applying  the  divergence  operator  to  Eq.  (12),  assuming  V x 0. 
The  result  in  normalized  units  is 

r+'ir  = V2[(r2  - l)cr]  + r + o-.  (13) 
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We  took  for  a solution  in  the  form  <r  — 5F (/)  W(r)  Yh„(b,  <£),  where  Ylm  is  the  usual  spherical 
harmonic.  Two  ordinary  differential  equations  result: 


'^T\7r  ((r2  - 1)  W)  - HI  + l)(r2  - 1)  W 
dW 

+ /■  — + (M  + 1)I^  = 0 
dr 


f"+23  = - aiJ, 


where  /u  is  a separation  constant. 


Equation  (14)  is  solved  together  with  the  boundary  conditions  fF(r±)  = 0.  First  we  show 
that  the  eigenvalue  m must  be  positive.  Letting  U = (r2  - 1)  W-'and  x = r 2,  we  can  rewrite  Eq. 
(14)  in  the  Sturm-Liouville  form 


4 -f  a3/j(x  - 1) 


v-i  dU  , u + 1 


I l y - 1 


2*l/20r  - 1) 


A i 2 

- Jr3/2(jf  - l)*"1  - /(/  + 1)|(/  = 0.  (16) 

Multiplying  Eq.  (16)  by  (/and  integrating  from  1 to  r\  in  the  usual  way  yields  an  expression 
for  n as  the  ratio  of  two  positive  quantities,  whence  ^ > 0. 

The  radial  eigenfunction  can  be  written  as  a superposition  of  two  independent  integrals  of 
Eq.  (14),  expressible  in  terms  of  hypergeometric  functions  2F ,.  For  example 


W(r)  - C+r'2Fl  \ (b  + + d+).  \ (b  + - d+ );  2 + — l— ; 1 - r2 

l I y — 1 


+ j/-,  ^ (6_  + rfj.  i (6.  - </J;  2 + — 1 - r2 


*±-2  + -iT  ± / + | . 
y - 1 2 
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d± 


> - 1 4 


/ + ± 

1 + \ T 4 + — - — 1 

2 

2 y-l) 

1/2 


(19) 


normalization  condition. 


evaluation  of  n in  general  must  be  carried 


tional  principle  associated  with  Eq.  (15).  Approximate  results  can  be  obtained  analytically  in  a 
number  of  interesting  limits,  however.  Thus  in  the  thin-shell  limit  the  WKB  approximation 


yields  /i  = 


n n 


rjr,  - 1 


2 

• , n - 1,  2,  ....  For  thick  shells  (r+  — <»)  the  ground  state 


eigenvalue  found  using  WKB  is  n =2(y  — 1 ) / (/  + 1).  The  spectrum  can  readily  be  shown  to 
have  the  following  properties  in  general:  (1)  It  is  discrete,  contains  an  infinite  number  of  modes 
and  is  unbounded  above;  (2)  the  lowest  level  (minimum)  value  of  n increases  with  increasing  / 
artd  y and  (3)  decreases  with  increasing  shell  thickness,  as  also  does  the  level  separation. 
Clearly  which  depends  on  the  value  of  the  indices  / and  n (but  not  m)  is  just  the  dimen- 
sionless characteristic  frequency. 


Given  two  independent  solutions  of  the  time-dependent  Eq.  05)  can  be  written 


§</> 


— (/"'  - 1) 
a 


1/2 


r i 3 ,2  + i£>  3 2 - m.  3.  . 

2ai  It  + ~ - t + — : — ; 1 ~ / 


4a 


4a 


(20) 

(21) 


where  A - [8Ma  - (a  + 2)2]1/J.  At  / - 0,  5(0)  - 1.  5(0)  - 0 and  §(0)  - 0,  §(0)  - 1.  As 
f — 0,  the  standard  asymptotic  formulas  yield 


50)  ~ 


r — r 

/A 

2 

2a 

r T + ' 7 A r 

4 4a 

1 2 + /A 

4 4a 

(22) 
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§(/> 


VF 


/A 

2a 


3 ! 2 - /A 

2 

2 + /A 

4 4a 

4 

4a 

■ »+2  + iA 

/ 4 


+ C.C. 


(23) 


Evidently,  the  relative  amplification  57./  ( or  §//)  consists  of  a rapidly  oscillating  part  multiplied 

j <«-2)  \ (y-5/3) 

by  / = / . If  y < 5/3  (y  > 5/3)  this  factor  diverges  (vanishes)  as  / — * 0.  For 

y = 5/3,  the  perturbations  and  the  size  of  the  shell  go  to  zero  at  the  same  rate.  The  constants 
multiplying  the  factors  containing  / depend  weakly  on  p.  through  A,  however,  both  decreasing 
monotonically  as  p increases.  Figure  1 displays  the  absolute  values  of  these  constants  as  func- 
tions of  m for  a = l(y  = 4/3).  Since  p increases  with  both  / and  «,  lower-order  modes  have 
larger  amplitudes  then  higher  ones,  hence  are  more  disruptive. 

We  note  that  by  Eqs.  (6)  and  (7),  p — /-3y,  and  so  the  limiting  dependence  of  Eqs.  (22) 
and  (23)  on  pit)  is  precisely  that  given  by  Eq.  (1).  Of  course,  to  study  the  temporal  develop- 
ment of  the  perturbations  in  detail  it  is  necessary  to  integrate  Eq.  (15)  [or  evaluate  Eqs.  (20) 
and  (21)]  numerically.  This  presents  no  difficulty  until  times  when  the  oscillation  period 
becomes  very  small,  at  which  stage  the  asymptotic  formulas  can  be  used. 

Case  (i) ' (y  = 1)  is  treated  in  very  similar  fashion.  Equation  (13)  is  replaced  by 


f1  ir  =>  VV  + r + a, 
dr 


03') 


which  again  separates.  The  space-  and  time-dependent  equations  are  both  solvable  in  terms  of 
confluent  hypergeometric  functions.  The  spectrum  obtained  from  the  former  is  qualitatively  like 
that  found  for  y > 1 in  all  important  respects.  The  solutions  of  the  time-dependent  equation, 

/25F  = -My.  05') 

may  be  written  as 
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*(j '•  t;i°4 


where  d>  and  ¥ are  the  standard  solutions  of  the  confluent  hypergeometric  equation.  As  / — 0, 
the  asymptotic  forms  of  the  latter  yield 


/ (In  /) 


r 7* 


K(f)  ~ On  f)  2 . (23  ) 

Thus  only  the  former  solution  gives  rise  to  an  amplification,  with  ^(r)// diverging  as  a power 

of  In  f 

Turning  next  to  example  (ii),  the  uniform  density  case,  we  begin  by  noting  that  for  an 

1 ( Sp o I ( dPo  I 

anisentropic  unperturbed  state,  the  terms  — V 1-^-  f ' Vj  - V • 1-^—  f Vs  must  be 
added  to  the  right  hand  side  of  Eq.  (12).  In  scaled  variables,  the  equation  for  ( is  now 

f"+2l  - V[(rJ  - 1)V  (]  + V£-  r - rV  ■ f.  (24) 

It  is  no  longer  possible  to  obtain  a solution  with  V x £ — 0.  Instead  we  find  two  coupled 
equations. 


/"+2cr  - VJ|(/-2  - 1) tr]  - 2cr  - r V x w 


/,,+J  w,—  — + r x Vtr. 
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Elimination  of  co  between  Eqs.  (25)  and  (26)  yields  an  equation  for  <t  which  is  fourth  order  in 
time  and  second  order  in  space: 


f"+2(J“+2a )"  = f"+2 


V2[(,-2  - l),yi  - 3 or 


+ v2[(r2  - l)<r]  - r2V2<r  + rr.VVcr  + 2r  • V<r  - 2tr.  (27) 

A solution  is  again  found  by  separation  of  variables,  with  the  time  dependence  given  by 
Eqs.  (15)  and  (15  ),  and  the  radial  dependence  expressible  in  the  form  Eq.  (17).  The  depen- 
dence on  (i  of  d2±  is  now  nonlinear.  The  solution  of  the  eigenvalue  problem  leads  to  two  values 
of  m,  one  positive  and  one  negative,  of  which  only  the  f.  ner  can  give  rise  to  amplification.  No 
special  treatment  is  required  in  the  limit  y — 1,  wlvch  is  now  nonsingular.  The  spectrum  is 
qualitatively  similar  to  that  found  for  case  (i).  The  only  changes  in  the  time  history  of  a pertur- 
bation are  those  ascribable  to  the  changed  values  of  /x.  It  is  clear  that  this  time  dependence  is  a 
general  property  of  the  homogeneous  self-similar  motion  on  which  the  present  examples  are 
based. 


Thus  in  the  two  examples  treated  in  detail  here,  sound  waves  standing  in  an  imploding 
spherical  shell  can  become  large  in  amplitude  compared  with  the  diminishing  shell  dimension, 
at  the  same  rate  in  each  instance.  From  this  it  may  be  concluded  that  the  result  is  insensitive  to 
the  form  of  the  density  profile.  Since  the  time  dependence  agrees  asymptotically  with  the 
model-independent  estimate  of  Eq.  (1),  it  is  unlikely  that  the  shape  of  the  pressure  pulse  plays 
an  important  role.  Because  the  long-wavelength  modes  have  the  fastest  relative  amplification, 
neither  viscosity  or  other  dissipative  processes  are  likely  to  provide  a way  of  controlling  them. 
Thermal  conduction,  which  reduces  the  effective  y,  only  makes  matters  worse.  Layering  might 
help  by  raising  the  natural  vibration  frequencies.  The  tendency  of  all  materials  to  behave 
incompressibly  (y  » 1)  at  high  densities,  when  the  Fermi  energy  becomes  large,  means  that 
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the  relative  amplification  cannot  really  diverge.  To  determine  for  a given  initial  vibration  level 
the  actual  limit  on  aspect  ratio  or  compression  imposed  by  the  mechanism  discussed  here 
requires  calculations  using  a more  detailed  model. 
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